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. ABSTRACT 

■ In this paper we present X-Shooter time resolved spectroscopy and RXTE PCA light curves of 

the ultra-compact X-ray binary candidat e 4U 06 14-H09 1 . The X -Shooter data are compared to 
the GMOS data analyzed previously bv iNelemans et al] (|2004|) . We confirm the presence of 
cm and O II emission features at ^ 4650 A and ^ 5000 A. The emission lines do not show 

C/3 , evident Doppler shifts that could be attributed to the motion of the donor star/hot spot around 

■ the center of mass of the binary. We note a weak periodic signal in the red-wing/blue-wing 

flux ratio of the emission feature at » 4650 A. The signal occurs at P = 30.23 + 0.03 min in 
the X-Shooter and at P - 30.468 + 0.006 min in the GMOS spectra when the source was in 

^ ■ the low/hard state. Due to aliasing effects the period in the GMOS and X-Shooter data could 

CsJ , well be the same. We deem it likely that the orbital period is thus close to 30 min, however, as 

■ several photometric periods have been reported for this source in the literature already, further 
OO \ confirmation of the 30 min period is warranted. We compare the surface area of the donor star 

and the disc of 4U 0614-1-091 with the surface area of the donor star and the disc in typical 
psj ■ hydrogen-rich low-mass X-ray binaries and the class of AM Canum Venaticorum stars and 

argue that the optical emission in 4U 0614+091 is likely dominated by the disc emission. 

■ Assuming that the accretion disc in 4U 0614+091 is optically thick and the opening angle 

of the disc is similar to the one derived for the hydrogen-rich low-mass X-ray binaries we 
conclude that the donor star may in fact stay in the disc shadow. Additionally, we search for 
periodic signals in all the publicly available RXTE PCA light curves of 4U 0614+091 which 

' could be associated with the orbital period of this source. A modulation at the orbital period 



H . with an amplitude of ^ 10% such as those that have been found in other ultra-compact X-ray 

^ ■ binaries (4U 0513-40, 4U 1820-30) is not present in 4U 0614+091. We propose that the 

X-ray modulation found in 4U 0513-40 and 4U 1820-30 could be a signature of an outflow 
launched during the high/soft state near the stream impact region. 
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1 INTRODUCTION 

The source 4U 0614+091 is a persistent ultra-compact X-ray bi- 
nary (UCXB) candidate that belongs to the class of low-mass X-ray 
binaries (LMXBs). UCXBs are characterized by very short orbital 
periods, typically less than 80 min. The consequence is a secondary 
star with such a high density that it excludes main-sequence type 
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donor star jNelson et alJl986l:ISavoniie et aljl98^ . The donor star 
in 4U 0614+091 is most probably a CO white dwarf, based on the 
emission lines of oxygen and carbon found in the optical spectra 
and the lack of hyd rogen and helium lines ( INelemans et al.l l2004l: 
IWemer et al1l2006l) . An additional hint towards a CO white dwarf 
comes from the de tection of the relativistically broadene d O VIII 
Lyff emission line jMadei et al1l2010l : ISchulz et al.ll2010t) . The ac- 
cretor is a neutron star, since it is showing type I X-ray bursts 
jKuulkerset a l. 20Q9J. The source shows no eclipses or dips in 
the light curves. However, due to a very low mass ratio (theoret- 
ical estimate based on the suggested orbital period of 50 min is 
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Figure 1. The average normalized X-Shooter spectrum of 4U 0514+091. Interstellar features, diffuse interstellar bands and telluric absorption lines are marked 
with 'i.s.', 'dib' and a horizontal line, respectively. Top panel: The average UVB arm spectrum. Note the emission feature at » 4650 A which most probably is 
a blend of carbon and oxygen lines. The features around 3800 A, marked with a dashed line are artifacts and they occur due to the imperfect merging of the 
master flats. Bottom panel: The average VIS arm spectrum. The residuals from the sky subtraction become more prominent above 7000 A. 



Table 1. Log of the X-Shooter observations. The time column indicates the start time of the first and the last observation obtained during each night. The third 
column shows the number of exposures in each arm. The average and ranges of the seeing values as well as the airmass during the time when the observations 
were taken are indicated. The seeing was measured from the spectra in the wavelength range 4500-4520 A. The range of S/N in the single spectra is calculated 
at a; 4500 A (UVB), a 6500 A (VIS). In the case of the NIR arm we provide the S/N of an average spectrum for each observing night, measured around 12500 

A. 



Date UTTime #ofexp.* Seeing ["] aver., range S/N UVB S/N VIS S/N (aver.) NIR Airmass 



16/02/2010 


01:33:32-02:20:01 


15 


0.81,0.73-0.96 


3.8-4.3 


5.2-6.3 


1.2 


1.21-1.26 


18/02/2010 


00:38:18-02:30:10 


30 


0.92, 0.66-1.12 


5.4-7.1 


6.9-9.8 


0.8 


1.21-1.29 


20/02/2010 


00:36:55-01:30:42 


15 


0.93,0.79-1.04 


4.0-7.2 


5.4-9.4 


1.3 


1.21-1.21 


21/02/2010 


00:29:50-01:23:37 


15 


1.03,0.91-1.21 


5.8-6.6 


7.4-8.7 


1.3 


1.21-1.21 


22/02/2010 


00:30:39-01:24:26 


15 


0.94, 0.85-1.05 


6.0-6.9 


7.7-8.9 


1.5 


1.20-1.22 



* In the case of the NIR arm three exposures were taken and averaged during the time of one UVB and VIS exposure 



q ~ 0.01) the threshold for the inclination below which the e clipses 
and dips are not prese nt can even be higher than 80 deg jHomd 
11985*: 'Paczvnsk illl974h . 

The orbital period measurement is important to classify the source 
as an UCXB and determine its evolutionary stage. In the case of 
4U 0614-1-091 th e first claim of t he orbital period measurement of 
Si 50 min was bv lO'BrienI j2005l) in high-speed optical data taken 
with ULTRACAM. The modulation in the optical light curve can 
occur either due to heating of one side of the donor by X-ray irradi- 
ation, the superhump phenomenon caused by a tidal deformation of 
accretion disc, emission from the impact point of the stream with 
the disc or an asymmetry of the accretion disc. The orbital period 
of ~ 50 min was also sugg ested by Shahbaz et al.l ( l2008h in IAC80 
data and Zhang et al. jlOlZ) in Otto Struve Telescope data. How- 



ever, iHakala et alj ( 120 Uh analyzed a large set of fast optical pho- 
tometry obtained with the NOT telescope and report no conclusive 
evidence for the orbital period of this system. In one out of 12 opti- 
cal data sets the authors found a strong periodic signal at » 50 min. 
The modulations with the orbital periods in the UC XBs were re- 
ported not only in optica l but also X-ray light curves jpiocchi et al.l 
l201lHStellaetalJll987h . As far as 4U 0614-1-091 is concerned a 
period ic modulat ion of x 25 min was found in one EXOSAT light 
curve faakala e "t al. 201 1) , when the sou rce was in a high flux state. 
On the other hand, Schulz et al. folded a Chandra HETGS 

light curve, when the source was in the low/hard state, with the pe- 
riods of 41, 51, 62 and 120 minutes and found no evidence of a 



periodic modulation. 

A powerful method of measuring the orbital parameters in persis- 
tent LMXBs or LMXBs in outburst is to measure the Doppler shifts 
of the high-excitation emission lines of C III and N III forming 
the Bowen blend at ^ 4650 A dSteeghs & Casaresll2002l) . These 
narrow lines are associated with the irradiated companion st ar and 
have been detected in memy LMXBs (Casaresetal. 2003. l2006l : 
IComelisse et al.ll2007l : iMurioz-Darias et al.ll2009l) . 
Another method used to determine the mass ratio is to measure the 
rotational broadening of th e absorption lines ori ginating in the at- 
mosphere of the donor star ( Gie s & Bolton|[l98^ . 
Nelemans et al. (2006) analyzed a set of time-resolved spectra 
of 4U 0614-)-091 taken with Gemini Multi-Object Spectrograph 
(GMOS) on the Gemini-North telescope and found no signatures 
of narrow emission lines. Although the source is showing a broad 
emission blend of O II and C III lines at a 4650 A the resolution of 
» 2 A might have been too low to detect narrow line components, 
if present. The authors found, however, weak evidence for a sinu- 
soidal pattern with a period of 48.547 min in a weak absorption line 
around 4960 A (possibly caused by C I), which could result from 
the absorption in the atmosphere of the donor star 
A spectroscopic measurement of the orbital period has not been 
successfu l for any of the UCXB s except for the weak evidence 
noted by iNelemans et alj ( |2006|) . It has been, however, success- 
ful for AM Canum Venaticorum (AM CVn) stars. This class of 
sources consists of a white dwarf accretor with a degenerate or 
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semi-degenerate, hydrogen-poor donor star. They are similar to 
UCXBs in terms of the mass ratio and chemical composition of the 
donor star. The Doppler shifts of emission lines found in the optical 
spectra of AM CVns are caused by the motion of the impact point 
of the accretion stream o nto the accretion disc (e.g. iNather et al.l 
ll98lH Roelofs et al. II2OO?) or, in the case of the ver y short period 
system HM Cnc, the irradiated side of the donor star jRoelofs et al.l 

In this paper we analyze time-resolved X-Shooter spectra of 
4U 0614+091 searching for narrow emission/absorption lines 
showing periodic Doppler shifts associated with the donor star/hot 
spot motion around the center of mass. The X-Shooter resolution of 
« 0.9 A in the blue part of the spectrum gives us a better opportu- 
nity to detect the narrow emissio n components t han previous stud- 
ies using the GMOS instrument iNelemans et al. 2006}. We com- 
pare the X-Shooter and the GMOS data. Additionally, we analyze 
all archival RXTE PCA light curves in search for an X-ray modu- 
lation at the orbital period. We separate high and low X-ray flux 
regimes of the source and look for periodic signals in both of them. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 X-Shooter observations 

4U 0614+091 was observed with the X-Shooter instrument on 5 
nights from 16-22 February 2010 (see Table. 1). X-Shooter is an 
echelle spectrograph operating at the Cassegrain focus of Unit Tele- 
scope 2 (Kueyen) on the Very Large Telescope. In a single exposure 
the instrument cov ers the spectral ran ge from the UV to the near-IR 
K band in 3 arms jVernet et al]|201lh . The ultraviolet-blue (UVB) 
arm covers 3000-5595 A in 12 orders with a resolution of 5100 (for 
a 1" slit width), the visual-red (VIS) arm covers the range 5595- 
10240 A in 14 orders with a resolution of 8800 (for a 0.9" slit 
width) and the near-infrared arm covers 10240-24800 A in 16 or- 
ders with a resolution of 5300 (for a 0.9" slit width). The slit widths 
described here were used during the observations. 
The source was observed in the slit stare mode (fixed position on 
the slit) in all three arms. The binning factor for the UVB and VIS 
arm spectra was 2 in the spatial direction and 2 in the spectral di- 
rection, the NIR arm spectra are by default not binned. The slow 
readout mode with high gain (100k/lpt/hg/2x2) was used for all 
the nights except the first night when the readout mode was fast 
with low gain (400k/lpt/lg/2x2). The exposure time of every ob- 
servation was 180 sec in UVB and VIS arm. Due to significant 
thermal background noise in the NIR arm instead of one 180 sec 
exposure, three exposures of 60 sec were taken and averaged in one 
FITS file. We reduce the data using the standard X-Shooter reduc- 
tion p ipeline version 1.5.0 ( Goldoni et al. 200^ : iModigliani et al.l 
I2OIOI) . A physical model mode was used during the reduction pro- 
cess in which the solution is obtained by optimizing the instrument 
physical model parameters to the data. Since the data were ob- 
tained in February 2010 during the first, shared-risk semester of 
X-Shooter operations, no linearity frames are available for these 
observations, therefore we do not compute the detector gain and 
linearity (xsh_lingain recipe). The bad pixel map is included only 
during the NIR arm data reduction process, since it contains signifi- 
cantly more bad pixels than the UVB and VIS arm data. In the UVB 
arm we use an instrumental response updated by the X-Shooter cal- 
ibration team which contains a more accurate than currently avail- 
able in the pipeline flux calibration around 3600-3800 A (Wolfram 
Freudling, private communication). In the VIS arm we use the in- 
strumental response file provided with the pipeline. The calibration 
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Figure 2. Fluxes averaged over one observing ID of all RXTE PCA, PCU 2 
observations in the X-ray energy ran ge 2-5 keV. Th e blue and orange mea- 
surements are taken from van Straate net alj )2000 j) and represent a sample 
of light curves in the island state up to lower banana state (blue) and lower 
banana up to higher banana state (orange). The red solid line at 70 cts/sec 
distinguishes roughly between the high/soft and low/hai'd state. 



of the source spectrum is additionally corrected for the instrumental 
flexures during the observing nights using the xsh_flexcomp recipe. 
In order to estimate the sky background, a one-dimensional spec- 
trum of the sky is built. The pixels outside of the source extraction 
region and far from the slit edges are taken. The interpolation of 
the resulting one-dimensional sky spectrum is done by applying 
the running median with the default window width of 14 pixels. 
The optimal extraction is used to obtain the one-dimensional spec- 
trum of the source. 

The analysis of the spectra is performed using the molly package 
developed by Tom Marsh. 



2.1.1 Average spectrum 

We extract all UVB and VIS spectra from every night and calculate 
the weighted (based on S/N ratio) average UVB and VIS spectrum 
(see Fig. 1). The resolution in the output UVB and VIS spectra is 
X 0.9 A. All the spectra are extracted with a wavelength sampling 
of 0.3 A and a median smoothing is used with a window width of 4 
bins for UVB arm and 5 bins for VIS arm. The plate scale is a 0.32 
arcsec/pix for UVB and VIS. We determine the seeing for every 
observation by measuring the FWHM of the source profile of each 
spectrum in the wavelength range 4500-4520 A (see Table 1 ; for the 
seeing measurement the spectra are unsmoothed and binned to the 
resolution of the instrument). The accuracy of the wavelength cali- 
bration in every spectrum is determined by measuring the position 
of the O I sky lines in all spectra. The sky line at 5577.338 A has 
an rms scatter of 0.05 A (UVB arm) and the sky line at 6300.304 A 
has an rms scatter of 0.07 A (VIS arm). 

Since one of the goals of this paper is to look for periodic Doppler 
shifts in the emission lines that could correspond to the orbital pe- 
riod of this source, we examine the spectra 'by eye' in search of 
prominent emission lines. We co nfirm the presence of an emi ssion 
feature at ^ 4650 A reported bv INelemans et alj ( l2004ll2006l) and 
IWemer et al.l (12006') corresponding to the blend of C III and O II 
lines (see Fig. 1). We also find weak emission line signatures of O II 
at a 4900 A and C II at 7200 A jNelemans et al.l2006l ; lw"emer et al.l 
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Table 2. The measured equivalent widths for the emission features: C III 
and O II transitions at 4624-4680 A and O II transitions at 4680-4720 A 
and 4900-4960 A. 



Date 


EW (C III, O II) 


EW (O II) 




(4624-4680 A) 


(4680-4720 A, 






4900-4960 A) 


16/02/2010 


6.7±0.1 


4.0±0.2 


18/02/2010 


3.25±0.06 


1.83±0.09 


20/02/2010 


4.67±0.09 


2.4±0.1 


21/02/2010 


5.59±0.1 


2.8±0.1 


22/02/2010 


4.6±0.1 


2.5±0.1 





0.4 



0.1 
0.0 



I I I 




I2OO6I ). The combination of two different flat fields (D2 and QTH) 
causes a discontinuity in the response curve and results in the arti- 
facts in the flux calibrated spectra around 3800 A hence this region 
is difficult to examine for possible emission features. Therefore, 
we chec k individual orders in search for emission lines of O III 
mostly ( Werner et al. I l200^ around 3700 A. However, we find no 
significant emission features that could be attributed to the emis- 
sion lines of oxygen or carbon in that region. The blue part of the 
UVB spectrum (3000-3600 A) shows small artifacts ('bumps' and 
'dips') visible also in the spectrum of the standard star possibly 
caused by imperfect merging of the spectral orders in the region 
where the S/N ratio of the data is much lower than in the middle 
and red part of the spectrum. 

We extract the median NIR spectrum from all observations taken 
each night. The resolution of the NIR spectrum is » 3 A and 
plate scale is » 0.21 arcsec/pix. We detect emission at wavelengths 
where the near-IR J, H and K bands occur. Most of the LMXBs are 
located towards the Galactic center. Therefore, their optical light is 
usually highly absorbed. A near-IR spectrum could prove an alter- 
native way to study the chemical composition of the donor star of 
the UCXBs and distinguishing between the CO or He WD. For that 
reason we look for significant differences between the LTE model 
spec tra for UCXBs with a carbon-oxygen and helium WD donor 
star jNelemans et ai]|2006l). There is a strong emission line at a; 
1.2 nm (blend of C III transitions: 1199.13, 1198.8 and 1198.12 
nm) present only in the CO model. Therefore, its presence could be 
used to distinguish an UCXB with a CO WD donor from He WD 
donor star. We examine the wavelength range where this emission 
line should occur. We find, however, no evidence for the presence 
of such a line. The upper limit on the flux of this line using the 
spectrum is also difficult to determine due to the systematic errors 
in the flux calibration. Due to a poor S/N ratio of the NIR spectra 
and many residuals from the sky lines subtraction we do not discuss 
this part of the spectrum further. 



2.2 RXTE observations 

We choose all of the Standard 2 mode (FS4a*.gz), PCU 2 data ob- 
tained during the /JXr£ mission (» 16years)of4U0614-)-091. The 
observations where type I X-ray bursts were detected (2 light curves 
which contain tails o f the X-ray bursts at MJD=51789 and 51944, 
iKuulkers et aLll2009i) are excluded and the ftools package is used 
to process the rest of the RXTE data. The data are screened with 
standard criteria and the light curves are extracted using the saex- 
TRCT tool in the channel range 2-14 (roughly 2-5 keV[J) with time 
bins of 16 sec. In this way our analysis is more sensitive to mech- 
anisms causing the modulation of the soft X-rays. The background 
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Figure 3. Top panel: RXTE ASM and MAX! flux of 4U 0614+091 close 
in time to the X-Shooter observations. The red line shows the approxi- 
mate threshold between the low/hard and high/soft state determined from 
the PCA fluxes (see Sec 3.4). The blue vertical lines indicate positions of 
the X-Shooter observations. A flux of 1 Crab=2.4xl0"'* erg/sec/cm^ corre- 
sponds to X 75 cts/sec/SSC for the ASM and 3.6 ph/sec/cm- for the MAXI 
instrument. Bottom panel:. RXTE ASM flux of 4U 0614-1-091 close to die 
optical CMOS observations (blue vertical lines). 



light curve is created by the rltnpcabackest tool using the back- 
ground model for bright sources and subtracted from the source 
light curve using the lcmath tool. We do not apply any correction 
for the deadtime. In total there are 488 RXTE PCA obsIDs and the 
length of the individual observations is less than 55 min. The major- 
ity of the observations have a length in the range 40-60 min which 
is similar to the orbital period proposed so far for 4U 0614-1-091. 
Hence if the orbital period of 4U 0614+091 is indeed a 50 min we 
may expect prominent aliases in the periodogram which can make 
the detection of the true signal more difficult. 
The light curves are analyzed using the Period Time-Series Analy- 
sis Package available in starlink namaka package. 



3 ANALYSIS AND RESULTS 
3.1 Optical emission 

3. LI The structure and variability of the emission features 

In order to investigate if there are any changes in the mean 
wavelength of the emission feature caused by e.g. changes in 
the ionization state between the observations, we cross-correlate 
the average spectrum from every observing night covering the 
emission feature at ~ 4650 A with the average from all observing 

' The gain of the PCA instrument used to drift slowly with time, which 
slightly changed the energy to channel assignment. 
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Figure 4. Upper panel: Trailed spectrogram of all of the UVB observations 
showing the same blend of C III and O II as in the lower panel around 4650 
A. Black stripes separate the data from each night. The width of the emis- 
sion feature is FWHM ~30A~ 1910 km/sec and is consistent with being 
constant during the observations. Lower panel: The average UVB spectrum 
around « 4650 A showing the emission feature which includes C III and 
O II transitions, most likely broadened by the motion in the accretion disc. 
The constant and four Gaussian lines fitted to the feature are oveiplotted 
(red line). The black and green vertical lines indicate O II and C III transi- 
tions, respectively, identified bv lNelemans et al] j2004h in their best-fit LTE 
model. The light blue, light red areas below the spectrum indicate the flux 
range of the blue and red wing of the feature we use in the method to search 
for periodic signals, respectively. 

nights using xcor in molly. For the data analysis of the emission 
features at ~ 4650 A and ~ 4690 A we use the errors based on 
the variance of the data in the wavelength range 4500-4600 A. 
The cross-correlation reveals that the emission feature is stationary 
between the observing nights within the Icr errorbar calculated by 
the XCOR function. The width of the feature also does not change 
significantly between the observing nights. 

We measure the equivalent widths of the emission features between 
4624 - 4680 A and 4680 - 4720 A together with 4900 - 4960 
A. We note a significant decrease in the equivalent widths of the 
emission features between the first and second observing night 
which seems to correlate with the increase of the seeing value (see 
Table 1 and 2). Since the seeing conditions were variable over the 
observing nights and the signal to noise ratio of the spectra and 
emission features low it is difficult to say whether the variations in 
the equivalent widths of the emission features are caused by the 
intrinsic variations of the source emission. 

In order to estimate the width and centroid of the emission features 
at » 4650 A and x 4690 A we first fit two Gaussian lines to the 
average spectrum from all observing nights which give a poor fit 
with reduced chi-squared value of ^^/v = 1476/360. Additional 
visual inspection of the features lead us to fit four Gaussian lines 
(see Fig. 4, lower panel). The reduced chi-squared value of the 
fit using four Gaussian lines is x'/'^ = 761/360. The centroids 
and widths of the three distinct components forming the feature 



at X 4650 A are: ^^eni = 4640.7 + 0.1 A, FWHM = 9.3 ± 0.3 A; 
/icen2 = 4649 ± 2.0 A, FWHM = 5 ± 4 A; ^.^ns = 4659 ± 2 A, 
FWHM = 20 ± 4 A. The centroid and width of the O II emission 
lines at 4680-4720 A are = 4694 ±2 A and FWHM = 30±5 A. 



3.1.2 X-ray flux during the optical observation 

We check the RXTE ASM and MAXI flux during the optical obser- 
vation of the source (see Fig. 3). The flux is measured in the energy 
range 1.5-12 keV by ASM and we average 10 observations. In the 
case of the MAXI instrument we use the one day averaged flux in 
the energy range 2-20 keV. The red line indicates the approximate 
threshold between the low/hard and high/soft state based on the 
flux of the source (see Sec 3.4 for details) calculated for the energy 
range 2-20 keV (assuming power-law slope of F = 2 and interstel- 
lar absorption of Nh = 3 x 10"' cm"^). The average X-ray flux 
level is around 10"' erg/cm^/sec (2-20 keV energy range) and indi- 
cates that 4U 0614+091 was in the low/hard state during the optical 
observations. 



3.2 Trailed spectra 

iNelemans et al. liooi) noticed possible variability in the weak ab- 
sorption line around 4960 A (possibly C I at 4959 A) with a period 
of 48.547 min. Therefore we extract all UVB spectra separately 
and plot all the spectra in the form of a trailed spectrogram. We 
use the wavelength sampling of 0.3 A and median smoothing with 
a window of 4 bins. The normalization of the spectra is done by 
fitting a third order spline to the continuum of each spectrum with 
the most prominent carbon and oxygen emission lines masked. 
Each spectrum is then divided by the function resulting from the 
fit. 

The expected radial velocity semi-amplitude of the primary is 
A"! = 5 X sin i km/sec and that of the secondary is K2 = 727 x sin i 
km/sec (assuming a 1.4 Mq mass neutron star, 0.01 Mg mass 
donor star and orbital period of 50 min). The maximum velocity 
broadening due to the integration time of a single exposure is 
2nK2 X Texp/Z'orb ~ 270 km/sec. Here, i is the inclination of the 
orbital plane with the plane of the sky. 

We examine the region around 4959 A carefully and phase all of 
the spectra with periods in the range 48-49 min and step of 0.001 
min. We fi nd, however, no signi ficant modulation such as that 
reported bv lNelemans et alj ( l2006h . 

A useful method to find the orbital period in time-resolved 
spectra is to look for the periodic chan ges in the flux ratio of the 
blue and red wing of the emission lines jNather et alj|l98l | ). It has 
been successfully applied to AM CVn sources (e.g.l Roelofs et al.l 
2005). In our case most emission features are not single transitions 
of e.g. helium as is the case in AM CVns but blends of oxygen 
and carbon transitions (see Fig. 4, lower panel). Therefore the 
periodic signals detected using this method could be caused by the 
Doppler shifts of the emission lines as well as the variable flux 
ratios between the carbon and oxygen emission lines forming the 
feature. 

We apply this method to the most prominent feature at » 4650 
A. We measure the ratio R between the blue and red wing of the 
feature R = Flux(4635 - 4650A)/Flux(4650 - 4665 A) for every 
observation and calculate the Lomb-Scargle periodogram. The 
periodogram shows a maximum power at P = 30.23 ± 0.03 min 
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Figure 5. Left panel: Lomb-Scai'gle periodograms using X-Shooter (al & al panels) and CMOS (bl & b2 panels) data calculated from the flux ratio of the 
blue and red wing of the emission feature at 4650 A measured in each individual spectrum. Note the peak at « 50 cycles/day a 30 min in the X-Shooter (al 
& a2 panels) as well as CMOS data (bl & b2 panels). The X-shooter periodogram is overplotted over the CMOS zoomed-in periodogram (orange curve, b2 
panel). The P2 probabilities for the highest peaks in the periodograms and two aliases. Right panel: The trailed spectrogram constructed by folding the data 
on the period found in the Lomb-Scargle periodograms: P = 30.23 ± 0.03 min for X-Shooter data (upper panel), P = 30.468 ± 0.006 min for GMOS. The 
spectrogram it plotted in velocity scale with one bin corresponding to 19.35 km/sec/px at 4650 A for X-Shooter and 60.43 km/sec/px at 4650 A for GMOS. 
Two periods are shown in the spectrograms. Note a weak periodic signal in the left, middle and right side of the feature. The average spectrum and average 
flux at « 4650 A is shown in the panel below and on the right side of the trailed spectrograms, respectively. 



(see Fig 5., left panel). We note that choosing a different binning of 
the data (e.g. 0.4 A) and a different window width for the median 
smoothing (e.g. 3 bins) causes the highest peak in the periodogram 
to shift to the position of the closest alias (e.g. P = 29.61 + 0.03 
min). Hence it is likely that the true error on the detected period is 
larger than the quoted value. 

The significance of the period is calculated using a Monte 
Carlo method with the number of permutations set to 1000. A 
Fisher randomization ( Linnell Nemec & Nemec 1985) is used 
to determine the significance of the peak and two significance 
estimates are given in the output. The first is PI and represents 
the probability that, given the frequency search parameters, no 
periodic component is present in the data with this period. The 



second significance P2 represents the probability that the period is 
not equal to the quoted value but is equal to some other value. In 
our calculation the PI » 0.0 (PI lies between 0.00 and 0.01 with 
95% confidence level) and P2 = 0.027 + 0.005, which indicates 
that this period may well not be the true period. 
Next, we create a trailed spectrogram with spectra folded on 
an orbital period of 30.23 min. For this calculation the spectral 
wavelength scale is converted into a velocity scale with Av = 19.35 
km/sec/px at 4650 A using the vein function in molly. The he- 
liocentric correction is additionally applied to all of the analyzed 
spectra. The trailed spectrogram reveals a very weak periodic 
pattern (see Fig. 5, right panel) on the left, middle and right side of 
the emission feature. The pattern on the blue side seems to be in 
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antiphase with respect to the pattern in the middle and on the red 
side. L ooking at the LTE mo del fitted to the emission feature (see 
Fig. 4, iNelemans et alj2004h the structure on the left side consists 
of the O II lines whereas the structure in the middle and right side 
is dominated by the C III lines. 

Additionally, we create trailed spectrograms between 28-32 min 
with the step of 0.001 min. Since there are 10.000 frames to 
examine we create a movie ordering the frames with increasing 
orbital period. We notice a repeatable pattern separated by 1 cy- 
cle/day which correspond to the aliases seen in the Lomb-Scargle 
periodogram. Owing to the weakness of the signal it is difficult 
to attribute the periodic signal found to the orbital period of the 
source. 

We investigate also the range of periods between 45-55 min by 
creating a number of trailed spectrograms with the step of 0.001 
min. None of the frames, however, show evidence for a periodic 
signal that could corre spond to the orbital p eriods suggested so far 
for 4U 06I4-F09I (e.g. lShahbaz et alfcOOsI) . 



3.3 Comparison between X-Shooter and GMOS data of 
4U 0614-H091 

We reanalyze time-reso lved spectra of 4U 061 4-1-091 obtained by 
the GMOS-North (see INelemans et al.ll2006h . In total there are 
52 spectra taken between 30/10/2003 and 22/11/2003. The spec- 
tra are normalized and the wavelength scale is converted into ve- 
locity scale with Av = 60.43 km/sec/px at 4650 A. We calculate 
the Lomb-Scargle periodogram from the flux ratio of the blue and 
red wing of the emission feature at ~ 4650 A in the same way as 
for the X-Shooter data. The periodogram shows maximum power 
at P = 30.468 ± 0.006 min (see Fig 5., left panel) with probabilities 
PI X and P2 = 0.69 ± O.OI. The detected period is close to the 
one detected in X-Shooter data but not the same. We calculate the 
P2 probabilities for the three highest peaks in the X-shooter and 
GMOS periodograms (see Fig 5., left panel). The P2 probability 
has a lower value for the alias on the left hand side of the highest 
peak in the X-shooter periodogram which indicates that the value 
of the period is more uncertain than indicated by the error on the 
highest peak. A high values of the P2 probability for all three peaks 
in the GMOS periodogram indicates that the detected period is un- 
certain as well. 

Additionally, we create a trailed spectrogram by folding GMOS 
spectra with the detected period (see Fig 5., left panel). The pat- 
tern in the trailed spectrogram seems to be similar to the one found 
in X-Shooter data. 



3.4 RXTE light curves: modulation of the X-ray emission 

First, we consider all available RXTE data. In order to normalize 
the light curves we divide each by its average flux. Since the light 
curves are not evenly spaced in time, we use the Lomb-Scargle 
method to search for periodic signals. We calculate a single Lomb- 
Scargle periodogram taking the 16 year long dataset. The peri- 
odogram is calculated in the frequency range 10""* - 3 x 10"^ Hz 
(6-167 min). We choose a frequency resolution of A/ = 10"' Hz 
based on the overall (16 years) time interval of the data. 
The Lomb-Scargle periodogram reveals no clear periodic signal in 
the considered frequency range. The highest peak occurs at « 84.9 
min, which does not correspond to any orbital period proposed so 
far. Folding the data with this period gives a periodic signal with an 
amplitude of around 0.5%. We note that this periodic signal is close 



to the orbital period of the RXTE satellite l lWen et alj|2006l) . Since 
the X-ray modulation found in the UCXB 4U 0513-40 was de- 
tected in the high/soft state of the source, we divide the light curves 
into a high flux group with PCA 2-5 keV count rate > 70 cts/sec 
which indicates a high/soft state and low flux group with PCA count 
rate < 70 cts/sec which indicates a low/hard state. This rough divi- 
sion between the island (low/hard state) and ban ana state (high/soft 
state) is done based on the results obtained by Ivan Straaten etaU 
1 I2OOOI) concerning the correlation between the spectral shape and 
the timing features (indicatin g the state of the source) and source 
flux. The dataset analyzed bv lvan Straaten et all ( l200(i) covers ob- 
servations performed during the years 1996-1998. In Fig. 2 we 
mark in blue those observations that are in the island state based 
on the spectral behavior and in orange those that are in the banana 
state. We calculate the Lomb-Scargle periodogram for the low and 
high flux data sets. There is again no clear evidence of a periodic 
signal in the light curves covering separately high/soft and low/hard 
state neither around » 50 min nor around 20 - 30 min. The highest 
peak in the periodogram occurs atx 81.5 min in the low/hard state 
and X 78 min in high/soft state. Folding the data with this period 
gives a periodic signal with an amplitude of around 0.5% and 1.3% 
for the low/hard and high/soft state, respectively. 

4 DISCUSSION 

We have analyzed X-Shooter data and reanalyzed GMOS data 
of 4U 0614+091 and find no clear si gnature of the orbital 
period around 5 min suggested by Ishahbaz et alj ( l2008h and 
INelemans et al.l ( l2006i) . We notice a weak periodic signature in 
the X-Shooter and GMOS data around 30 min, which could well 
be due to the orbital period. On the basis of the strong aliases 
present around the most probable peaks in the X-shooter and 
GMOS periodograms we cannot exclude that the two periods are 
consistent with being the same. It is possible that the modulation is 
caused by the Doppler shifts of the C III and O II lines as well as 
variable flux ratio of the C III with respect to the O II lines forming 
the feature. Perhaps we are observing two regions characterized by 
different ionization states (e.g. disc and a stream impact region). 
In order to confirm this interpretation, however, spectroscopic data 
with a resolution higher than X-shooter in combination with higher 
effective area would be necessary in order to resolve the narrow 
components and obtain sufficient signal to noise ratio. 
Although there have b een claims of the orbital period ~ 50 
min using optical data jShahbaz etal] l2008l : IZhang et al] |20I2|) 
the authors found also indication for other periodic signals. 
IShahbaz et al.l 1 I2OO8I) reports two periodicities at 64.1 min and 42 
min in the NOT and SPM light curves, respectively. IZhang et al.l 
( I2OI2I) finds a possibly quasi-periodic signal also at I6.2±0.I min 
using data obtained by the Otto Struve Telescope. Therefore further 
confirmation of the signal at a; 30 min is warranted given the large 
number of reported periods and quasi-periodicities for this source. 

The lack of helium lines in optical spectra suggests that the 
donor star in the source is either a CO white dwarf or He star 
which is very late in its evo lution, on the (sem i-)degenerate branch 
of the evolutionary track dYungelsonl I2OO8I) . However, consider- 
ing the bolometric flux of 4U 0614-1-091 estimated from the RXTE 
ASM data the source appe ars too bright to have a degenerate donor 
jvan Haaften et 311120121/). The flux would match the theoretical 
pred ictions only if the orbital period of this source was around 25 
min jvan Haaften et alj2012bl) , which is more in line with the value 
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Figure 6. The theoretical donor star mass as a function of the orbital pe- 
riod for He (black) and CO WD (red). The solid curve represents a zero- 
temperature track which neglects the thermal pressure in the WD, whereas 
the dashed curve illustrates an upper limit on the donor mass in the case a 
very high thermal pressure is introduced. The dots indicate the logarithm of 
the age calculated since the onset of mass transfer to the neutron star, 

detected in the X-shooter and GMOS data than the claimed 50 min 
period. 



4.1 X-ray reprocessing: 4U 0614-1-091 compared to LMXBs 
and AM CVns 

Since we find no clear evidence for Doppler shifts in the velocities 
of emission lines produced on the X-ray heated hemisphere of the 
donor star we consider the possibility that the optical emission, in- 
cluding the emission lines, is dominated by reprocessing of X-rays 
and the contribution from the donor star to the total as well as the 
reprocessed emission is negligible. 

At first we estimate the visible area of the donor star and the disc 
as a function of the inclination. In this simple approach we assume 
that the ionization state and vertical structure of the accretion disc 
in hydrogen-rich LMXB and UCXB is similar. The shadowing of 
the donor by the disc is not included in this calculation, however, 
we discuss this problem later in this section. We use the formula 
describing the rat io between the do nor and the accretor Roche lobe 
radius derived bv lEggletonl ( ll983h . We estimate the projected area 
as a function of the inclination of the system assuming that the disc 
fills the entire Roche lobe of the primary 

5j ^ -^WD ^ Q.5nRl,jj(l + sinQ 

5disc + 5' WD ^^disc + 0.57ri?^j,(l + sin i) 

Sfrac is the ratio between the projected area of the WD (5wd) and 
the sum of 5wd and the projected area of the disc (5 disc)- We take 
into account in this calculation only the side of the donor star 
which is illuminated by the X-rays (e.g. for the inclination of 90 
deg half of the donor star total surface is illuminated and visible for 
the observer). The surface area of the donor star and the disc to first 
approximation depends on the mass ratio between the donor and 
accretor. In the case of the neutron star we assume the canonical 
mass of M = IAMq. We calculate the mass o f CO and He WD 
donor stars as a function of the orbital period jvan Haaften et al.l 
l20I2al lbl) assuming that the thermal pressure in the WD is zero 



Table 3. A list of LMXB and AM CVn sources for which the mass ratio 
was measured from spectroscopy. The inclination (if constrained) and the 
reference for each source is also provided. In the case of AM CVn we list 
only the sources with mass ratio close to that of 4U 0614+091, 



Name 


mass ratio 


inc, [deg] 


Reference 


LMXB 


ScoX-1 


0,3 


38 


Steeshs & Casares (2002) 


X 1822-371 


0,296-0,315 


82,5 ± 1,5 


Casares et al, (20(H) 


4U 1636-536 


0,21-0,34 


36-60 


Casares et al, (2006) 


4U 1735-444 


0,05-0,41 


27-60 


Casares et al, (2006) 


EXO 0748-676 


0,11-0,28 


75-82 


Muiioz-Darias et al, (2009) 


Cyg X-2 


0,34±0,02 


62,5±4,0 


Casares et al. (2010) 


LMC X-2 


<0,1 


<65 


Cornelisse et al, (2007) 


GX9+9 


0,07-0,35 


<65 


Cornelisse et al, (2007) 


AM CVn 


GPCom 


0,018 




Marsh (1999) 


SDSS J1240-01 


0,039 


36,53 


Roelofs et al, (2005) 


V396 Hya 


0,0125 




Ruiz et al, (2001) 



and the only forces present are the degeneracy pressure and the 
Coulomb interaction (Fig. 6, zero-temperature donor track). The 
Coulomb interaction becomes important in very low-mass WDs 
(< 10"^^ Mq). In order to give an idea what would be the donor 
star mass in the case where the thermal pressure is non-negligible 
we c alculate the donor star mas s with a radius larger by a factor 
of 2. iDelove & Bildsterj ( l2003h estimated that the influence of 
the thermal pressure on the radius of the white dwarf is probably 
around few tens of a percent. Hence, we should expect the donor 
star mass to be in between the solid and dashed curve (Fig. 6) and 
likely closer to the solid curve. 

Fig. 7 shows the resulting surface area of the white dwarf sec- 
ondary divided by the total surface area (donor star and the disc) 
as a function of the inclination (the mass ratio in 4U 0614H-091 is 
assumed to be g = O.OI). If the inclination of 4U 0614-1-091 is less 
than around 60 deg then the surface area of the donor with respect 
to the total surface area is significantly smaller for this source 
than for other, hydrogen-rich, LMXBs which show emission lines 
formed in the irradiated part of the donor star (see Table 3). Hence, 
the visible area of the donor and the disc in an UCXB with a low 
mass donor star and low inclination (less than around 60 deg) 
argues in favor of disc-dominated optical emission. 
The majority of t he UCXBs found so far have orbital periods 
around 40-60 min jin't Zand et ^l2007h . similar to that suggested 
for 4U 0614-1-091. The inclination of many is not well constrained. 
For most of them an upper limit is only given based on the lack of 
dips or eclipses in the light curve. However, since UCXBs likely 
have a lower mass ratio than hydrogen-rich LMXBs the inclination 
can reach val ues even higher than 80 deg without the occurrence 
of an eclipse jPaczynskUl 1 9741 : [Homell 1 9851) . On the other hand the 
probability of having systems with a very high inclination is quite 
low (see Fig. 7, top axis). 

Therefore, considering the difference between the visible fraction 
of the donor star and the disc surface in these UCXBs with respect 
to the LMXBs we may suspect the optical emission of UCXBs 
similar to 4U 0614-f091 to be dominated by the disc emission. 
As noticed earlier, the above calculation does not take into account 
the possible 'shadowing effect' of the disc on the donor star. 

© - RAS, MNRAS 000,[l]-?? 



X-Shooter spectra and RXTE light curves of the UCXB candidate 4U 0614+091 9 



P(InclInatIon>i) 

1.00 0.98 0.94 0.87 0.77 0.64 0.50 0.34 0.17 0.00 




20 40 60 80 

i [deg] 



Figure 7. The surface area of the donor star divided by the total surface area (donor star+accretion disc) as a function of binary inclination (see Eq. 1), 
oveiplotted are the best-estimate values for a number of LMXBs. Note that 4U 0614+091 could have a significantly lower donor star surface area compared 
to the total surface area than hydrogen-rich LMXBs for which the reprocessing of X-ray light in the donor star has been observed (red line). We plot (in blue) 
the positions of those AM CVns that have a mass ratio close to 4U 0614+091 and spectroscopic measurement of the orbital period. In the case we do not have 
the measurement of the inclination the source is plotted at 60 deg. The top axis represents the probability that the inclination of the system is higher than given 
for a homogeneous distribution in cos i. 



Hence, we estimate the height of the outer part of the disc and 
compare it with the theoretical estimates of the donor star radius. 
For a typical LMXB a half-opening angle of the disc is a 12 
deg dde Jong et al]|l996l) . The carbon-oxygen disc could, however, 
be thinner by a fact or of ^ 2 than a disc with solar composition 
in a typical LMXB ( IDunkel et al. 2006. whose prediction is based 
on the a-disc theory). The height of the outer part of the disc 
for a disc radius of 0.45 Rq is therefore in the range 0.045-0.09 
Rq. The radius of the donor star assuming the zero-temperature 
track is » 0.04 Rq which is still smaller than the lower limit on 
the height of the outer part of the disc. In the case in which the 
radius is twice the radius expected for the zero-temperature track 
(very high thermal pressure) it is possible that the donor is being 
illuminated. Such a case is, however, unlikely given the typical 
influence of the thermal p ressure on the radius of the white dwarf 
jPelove & BildsteiJl2003l) . 

Apart from the geometry of the accretion disc we consider also its 
optical thickness. The donor star could still be illuminated by the 
X-ray s that are passing through the accretion disc. Dunkel e t al.| 
( l2006l) have shown, however, that in the case of accretion disc 
dominated by ions of heavier elements such as helium, carbon or 
oxygen the electron and ion density as well as optical thickness of 
the disc increases (by a factor of few) with respect to the optical 
thickness of the disc in a typical hydrogen-rich LMXB. 
Based on the theoretical predictions presented above it seems 
likely that the donor star is shielded from the radiation of the 



central source. If that is the case then any modulation of the optical 
light is unlikely to be caused by the reprocessing of the X-ray light 
in the donor star as is the case for typical hydrogen-rich LMXBs. 
Fig. 7 indicates as well the position of a few AM CVn source for 
which the orbital periods were measured and the mass ratios are 
close to the predicted value in 4U 0614+091. Unlike in AM CVns, 
the accretion discs in LMXBs are strongly X-ray irradiated and 
the amount of reprocessing usually determines the temperature 
and elements visible in the optical spectra. Although the surface 
area of the donor star comparing to the accretion disc is similar in 
AM CVn source and UCXBs we may expect that the conditions 
in which the lines are formed in 4U 0614+091 are substantially 
different from those seen in AM CVn sources. 



4.2 X-ray modulation 

The origin of the X-ray modulation in the high/soft state that was 
observed for instance in the UCXB 4U 0513-40 ^ Fiocchi et al.l 
201 1) or UCXB 4U 1820-30 ( Stella et al. 1987) has not yet been 
established. In the case of the 4U 0513-40 a periodi c modulation 
is vis ible when the source is in the high/soft state jFiocchi et al.l 
l201lh . In 4U 1820-30 a periodic modulation with an amplitude of 
w 10% was found when the source was in the high/soft state. The 
modulation was also observed when the source was in the low/hard 
state, however, the significance of the detection is reduced due to 
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the factor of 3 lower count r ate and factor of 2 lower modulation 
amplitude jStella et al.lll987p. According to the model of repro- 
cessing presented bv lArons&Kinsl ( ll993h which describes the ob- 
servations of UCXB 4U 1820-30 the X-ray light is reprocessed 
into the optical and ultra-violet light mostly in the outer regions 
of the optically thick, geometrically thin accretion disc. Additional 
5-10% of the light comes from reprocessing by the degenerate 
donor and is causing a prominent modulation with orbital period 
of the optical and ultra-violet light in 4U 1820-30. However, the 
X-ray light reflected off the surface of the donor star cause a mod- 
ulation with an amplitude of only a few tenths of a percent of the 
total X-ray flux, hence it is insufficient to explain the large observed 
amplitude of the X-ray modulation in 4U 1820-30. 
We propose that the modulation found in 4U 1820-30 and 
4U 0513-40 may indicate the presence of an outflow launched 
during the high/soft state. The winds in several LMXBs are 
azimuthally sym metric and have an equatorial geometry (see 



uthatty 

IPonti et alj|2012f) . However, in the case of an UCXB the outer ra- 
dius to launch a thermal wind is comparable to the outer disk radius 
(a 10' cm. Miller et al. 2006). Therefore, when the point of impact 
of the gas stream from the companion occurs in the region where 
the wind is generated, the azimuthal symmetry is violated and more 
wind comes olf this impact region rather than off other regions in 
the disk. Hence the foreground absorption will depend on the or- 
bital phase. Additionally, since the wind in an UCXB is hydrogen 
deficient and helium or carbon and oxygen rich, the amount of ab- 
sorption may be higher due to the strong dependence of the photo- 
electric absorption cross-section upon the atomic number. 
A possible outflow p resent in 4U 1820-30 in the high/soft state 
has been reported by ICosta ntini et alj ( l2012h . The authors found 
highly blue-shifted absorption lines of mildly ionized oxygen. The 
observed lines could originate in the stream impac t region which 
could be less ionized than the accretion disc (e.g. I van Peet et al.l 
I2OO9I) . 

In order to study a dependence of the neutral and ionized absorp- 
tion on the orbital period using current grating instruments a long 
exposure time (wdays) is required. Taking into account the fact that 
the observed variability of the UCXBs occurs on a time scale of 
days this study will be more feasible using future higher-effective 
area instruments like micro-calorimeter on Astro-H satellite. 
Lack of a convincing periodic signal in the X-ray light curve with 
amplitude of few percent in the high/soft state of 4U 0614+091 
indicates that the mechanism causing the modulation of X-rays in 
4U 0513-40 is not present in 4U 0614+091. On the other hand the 
RXTE/PCA energy range, which does not cover the soft part of the 
spectrum (< 2 keV), could be the reason why the modulation re- 
lated to photoelectric absorption (most prominent in the soft part of 
the spectrum) is not detected. Alternatively a lower inclination of 
4U 0614+091 with respect to the UCXBs showing the X-ray mod- 
ulation could al so explain the lack of the periodic signal with such a 
high amplitude. IPonti et alj ( l2012h demonstrate that the winds ob- 
served in a sample of black hole X-ray binaries in high/soft state 
have opening angles of few tens of degrees since they are only ob- 
served in dipping sources in which the disc is inclined at a large 
angle to the line of sight (/ g 60 deg). A high-resolution LETGS 
spectrum of 4U 0614+091 taken during the high/soft state does not 
show signatures (na rrow absorption/emission lines) of an outflow 
( |Paerelsetal]|200lh . Therefore, assuming that the UCXBs share 
the outflow pro perties with the sa mple of black hole X-ray bina- 
ries analyzed bv lPonti et aD(|2012') the inclination of 4U 0614+091 
could in fact be lower than the threshold determined by the lack of 
eclipses in the X-ray light curve (see Sec. 4.1). 



The limit on the amplitude of any X-ray modulation in the RXTE 
data of 4U 0614+091 is such that reflection of X-ray light olf the 
surface of the donor star could still be below the detection thresh- 
old. If we assume that t he star atmosphere reflects a fraction w 0.5 
of the incident X-rays ( lAndersoij[l98ll : [London et al.lll98lh then 
the amplitude of the X-ray modulation is expected to be lower than 
1% (Arons & King 1993) for the case of an UCXB. However, the 
theoretical prediction seems to indicate that the disc is likely block- 
ing the X-rays from illuminating the donor star. 



5 CONCLUSIONS 

We have analyzed time-resolved X-shooter and GMOS spectra of 
4U 0614+091. We find a weak periodic signal at a 30 min in the 
red-wing/blue-wing flux ratio of the most prominent emission fea- 
ture at » 4650 A. The modulation could be caused by the Doppler 
shifts of the C III and O II lines as well as variable flux ratio of the 
C III with respect to the O II lines forming the feature. Comparing 
the surface area of the donor and the disc of 4U 0614+091 with 
the surface are of the donor star and the disc in typical hydrogen- 
rich LMXBs and AM CVn sources, we conclude that the emission 
likely originates in the accretion disc. It is possible that this peri- 
odic signal represents the orbital period of the source, however, due 
to large number of periods reported for this source a further confir- 
mation is needed. 

Additionally we find no evident periodic signal which could be at- 
tributed to the orbital period of 4U 0614+091 in the RXTE PCA 
light curves covering the energy range 2-5 keV. 
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